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NATIONAL ADVISORY COMITrEE FOR AERONAUTICS 

RESEARCH IvIEMORANDlM 

A SIMPLIFIED METHOD FOR APPROXIMATING THE TRANSIENT 
MOTION IN ANGLES OF ATTACK AND SIDESLIP DURING 
A CONSTANT ROLLING MANEUVER 
By Leonard Sternfield 

■ SUMMARY 


The transient motion in angles of attack and sideslip during a con- 
stant rolling maneuver has been analyzed. Simplified expressions are 
presented for -the- determination of the pertinenl^ modes of motion as well 
as the modal coefficient corresponding to each mode. 

Calciilations made with and without the derivatives Cyp (side force 
due to sideslip) and Cj^ (lift-ciarve slope) indicate that although these 
derivatives increase the total damping of the system, they do not markedly 
affect the transient motion. 


INTRODUCTION 


Recent flight tests of airplanes designed with their mass concen- 
trated primarily in the fuselage have indicated that dxuring a rolling 
maneuver the airplane experiences large uncontrollable motions in angles 
of attack and sideslip (refs. 1 and 2) . A fundmental analysis of this 
pitch-yaw divergence problem, assxmiing the rolling velocity of the air- 
plane is constant, is presented in reference 5- This analysis, concerned 
primarily with the stability of the system, makes possible the calciila- 
tion of the divergence boundaries and the prediction of the range of 
rolling velocities for which the airplane motions will diverge. Analog 
studies of several research airplanes reported in references A- and 5 
indicate that reference ^ Is a helpful guide in a roll-coupling investi- 
gation but detailed motion studies based on five degrees of freedom and 
taking account of pilot Inputs are essential to the analysis. Also, cal- 
culated time histories of the airplane motion indicate that large angles 
of attack and sideslip that are objectionable to a pilot and which may 
induce severe loads on the airplane may be encountered for constant 
rolling velocities outside the critical range. 
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The ptirpose of this paper is to extend the analysis of reference 3 
by analyzing the transient motion in angles of attack and sideslip. 
Simplified expressions are presented which permit an accurate and rapid 
estimate of the maximum angles of attack and sideslip for an airplane 
rolling at a constant velocity which shovild apply to the case of an air- 
plane entering a rolling maneuver up to the point of recovery from the 
maneuver. The airplane is assumed to be initially disturbed by an input 
in the side-force equation equal to the product of the rolling velocity 
and trim angle of attack. A comparison is made of the time histories 
obtained by using the simplified and exact expressions (obtained from 
the foirr degrees of freedom assuming constant rolling velocity) for a 
currently designed swept-wing fighter flying at a Mach number of 0.7 
at an altitude of 32,000 feet. 


SYMBOLS 

b wing span, ft 

c wing mean aerodynamic chord, ft 

Lif'b 

Ct lift coefficient, . 

^ ^ qS 

Cl rolling-moment coefficient, Rolling moment 

qSb 

Cm pitching-moment coefficient, Rllohing moment 

qSc 

Cn yawing-moment coefficient, 

qSb 

Ix moment of inertia of airplane about principal X-axis, slug-ft^ 

P 

lY moment of inertia of airplane about principal Y-axis, slug-ft 

I2, moment of inertia of airplane about principal Z-axis, slug-ft^ 

L' rolling moment, ft- lb 

m airplane mass, slugs 
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P 


PO 

q 

r 

S 

t 


'^1/2 


period, sec 

steady rolling velocity, radians/sec 

pitching velocity, radians/sec; dynamic pressure, Ib/sq ft 
yawing velocity, radians/sec 
wing area, sq ft 
time, sec 

time to damp to one-half amplitude, sec 


T 2 time to double amplitude, sec 

V airplane velocity, ft/sec 

Y side force, lb --------- _ . . . ' 

Oq initial angle of attack of airplane principal axis, radians 

Act. incremental change in angle of attack, radians 

a = oq + da 


|3 angle of sideslip, radians 

a + io) complex roots of characteristic equation 
At real root of characteristic equation 






lyPo^ 


BCni 

p — 
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, _ ^Cn 

SCl 


CLa = 




‘P ^ 




acj_ 

3b 
2V 

^Cm 

2V 




n 


' 2V 


^Oy 

> = i? 

La = <lSCx^ 


L'p = ^ C^ 

P 2V ‘P 


Ma = qScCni(j_ 
.2 


Mq = 


_ qSc 


2V 


^mq 


Np = qSbCnp 


Nr = 


qSb2 


2V 

Yp = gSCyp 


Cnj. 
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EQUATIONS OF MOTION 


In this analysis, assuming constant forward velocity and constant 
rolling velocity, the following linearized equations of motion of an 
airplane referred to the principal axes were used: 

Pitching : 


Yawing: 


Y-force : 


Z- force : 


q - 


~ ^ x' 
lY 


Mqq M^ 

PO^ = + — ^ 


(la) 


r ^ 


Iz 


Nv 


Nn 


POq = ^ r + -^ p- 

iz Iz 


(lb) 


P + r - po Al = _ p + PqOq 


mV 


(Ic) 


• . a ^ A 

o- - q + PnP ^ 

mV 


(Id) 


By setting q = q^e^^, r = rQe^^, p = PqB^^, and Ax = (Zkx)Qe^^ 

in the equations written in determinant form, A must be a root of the 
stability equation .... 


AA^ + + CA^ + DA + E = 0 


( 2 ) 
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where 


A = 1 


B = 


Nr _ ^ 

I^ ly mV mV 


C = 


ix - iy\/iz - ix\ p 


% % ^z 



iriV I^i iiiV ly 


L(x Nr Lg, Mg Lg, 

mV I^ mV ly mV mV 

/Nr ^ M,^ ^ MqNp^Mo^Nj. Mg N^ Yp ^ Mg N^ Lg^ ^ Np Lg^ ^ M^J^^ Yp ^ 

i Iz lyl ^ ly Iz ^Y .^Z ^Y ^Z ^Y ^Z ^Z ^Y 


i2 ,Y^%^is.Ya!!a + ^b Ax - IyV^Z - Ix' l 2 . Efi/Viixyjz^'lpo 

mV mV I'l mV mV ly mV^ Iz J\ ly / ^Z / \ ^Y / 


and 


E = - 


ix - I- 


xz -Ji^k - - ll\ 


Ix 


Po 


2 + 


Nr ^ Mq, Np Mq Np Jjq^ M(x Nr Yp Mq Nr La Yp ^ 

ly I2, ^ ly ^Y ^Y ^Z ^ ^Y ^Z 

LaY^&LlixVEL^ 

mV mV \ Iz / y ly / 
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In reference 5, Clq^ and were assumed to be equal to zero 

which corresponds to the case of the center of gravity moving in a 
straight line. In this analysis the case of CLq, = Cy^ = 0 is considered 

(referred to as case (a)) as well as the case where the values of Cx^ 


and 



are finite (referred to as 


case (b)). 


ANALYSIS AND DISCUSSION 


Divergence Boimdaries 


The conditions necessary for the system described by eqiiations (la) 
to (id) to be stable are that the coefficients.. of the characteristic equa- 
tion (eq. 2) must be positive and that Routh's discriminant, BCD - AD2 - 
B^E, be greater than zero. It can be shown from an examination of the 
characteristic equation of this system, that, for aircraft designed with 
positive-astatic -stability -and where each one of the damping-derivat ives - 
Cnr> Cniq, Cj^, and Cy^ adds damping to the system, the oscillatory 

modes will be stable for all values of the steady rolling velocity. Thus, 
the only type of instability that could be encountered is an aperiodic 
mode which occurs when the constant term of the characteristic equation, 
the E coefficient, is negative. As shown in reference divergence 
bomdaries, obtained by setting E = 0, can be plotted as a function of 
the squares of the natural frequencies in pitch and yaw nondimens ionali zed 
to the square of the rolling velocity, and respectively. 

These divergence boundaries are shown in figvire 1 for cases (a) and (b) 
for an airplane whose mass and aerodynamic characteristics are presented 
in table I. The nondimens ional pitch and yaw frequencies of the airplane 
described in table I fall along the straight line, with the slope eqml 

to — % shown in figure 1. Each point on this line corresponds 

Cnp ly b 

to a particular value of pq and as pQ increases the point moves along 

the line toward the origin. For case (a), the airplane line intersects 
the divergence boundaries at two points and thus defines the critical 
range of rolling velocities where a divergent mode exists, namely, 

1.86 < |pq| < 2.53. For case (b) the airplane line lies in the stable 
region and the motion is therefore stable for all values of Pq- 
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Characteristic Modes 


The roots of the characteristic stability equation corresponding to 
the values of jpQj shown in figijre 1 are presented in tables II (a) 

and II (b) for cases (a) and (b), respectively. A method of approximating 
the roots of the characteristic equation is presented in the appendix. 

The roots for pq = 0 which correspond to the condition where no coupling 

exists between the lateral and longitudinal modes are also given in the 
•table. If the roots are complex, A = a ± ica, P = — , and T-j_y'2 ~ ■ - ' * ' 


whereas, for a real root 


Ti/2 = 


-0.695 


^1 


or 


T2 


0.695 

^1 


The roots for pQ = 0 , for both cases (a) and (b), indicate that 
the motion consists of two independent oscillatory modes, a longitudinal 
and lateral oscillation with periods of 2."jk seconds and 4.08 seconds, 
respectively. The damping of the oscillation for case (b) is greater 
than that for case (a) because the damping derivatives Oy^ and CLq^ 

were assumed equal to zero in case (a). The values of Tj^/p in' seconds 
are given in the following table: 


Case 

^1/2 

Longitudinal oscillation 

Lateral oscillation 


(a) 

5.5 

15.2 

(b) 

1.4 

9.5 


For values of |pq| other than zero, the lateral variables p 

and r are coupled with the longitudinal variables a and q and the 
system now has the characteristic modes shown in table II for each one 
of the variables, p, r, a, and q. It should be noted that as |pq| 

increases, the frequency of the high frequency mode continues to increase 
and the damping decreases. However, the frequency of the other oscilla- 
tion decreases until |pq| reaches approximately 1.86' when the oscilla- 
tion breaks down into two aperiodic modes. In case (a) at this value 
of I Po I > S’ zero root is obtained and the unstable region begins . The 
system remains unstable, with one aperiodic divergent mode, until jpQj 

becomes equal to 2.55. In case (b), the system is stable although two 
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aperiodic modes occiir for 1.86 < |pq| < 2.53* For |pq| > 2.33^ "the 

oscillation reappears in cases (a) and (b) and the frequency and damping 
increase as |po| increases. 

From the contour lines of constant oscillation frequencies presented 
in reference 5, it appears that, in general, the frequency of one of the 
oscillatory modes will always decrease and probably break down into two 
aperiodic modes as IPqI increases. 


Motions in ^(t) and Zia(t) 


Time histories were calculated by the method of Laplace transform 
(ref. 6) of the motions in p and /to for the several values of Pq 

shown in figvire 1. The disturbance acting on the airplane is represented 
by the term Pqoo in the side-force, equation while the airplane is rolling 


in a negative direction. 


The motions in 


P(t) 

■°o 


and 


/to(t) 
‘ OQ- 


are presented 


in figures 2 and 5 for .cases (a) and (b), respectively. A comparison of 
figures 2 and 3 shows that the inclusion of and Cy^ modifies the 

motions slightly. Additional time histories were n^de, and are presented 
in figure for the same airplane but with the following assumed values 
of the derivatives: Cm^ = -0.09, Cn^ = 0.1l4, Cy^ = Cl^ = 0. For 

Pq = - 1.0, the airplane would be located in figure 1 at the point 

0)02 = 1*32 and o)^|,2 = 


Although the modes of motion Involved in the transient behavior of 
the airplane will be stable for points along the airplane line located in 
the stable region, that is, 2.33 < |Poj < 1-86 in case (a) and for all 

values of |po| in case (b), the magnitudes of /to and p developed 

during the transient motion may be large, enovigh to cause severe loads 
on the airplane. For example, in the particular flight condition for 
which the motions were calculated, where oq = 5*^^ airplane develops 

a sideslip angle of about -8^ at Pq = “1*5 end to. is about -20 at 
Pq = -3.0. These values are reached in about 3 seconds. 

In all the motions the high-frequency mode is clearly evident, 
although the amplitude of this mode is relatively small. Hence, the high- 
frequency mode modifies the shape of the motion only slightly. Thus, 
for Pq = - 1.0, -1.5, -2.5, and -3-0, the low-frequency mode is. of primary 

importance in the motions, whereas for po = -1*86, -2.0, and r2.33^ "the 
aperiodic modes are of primary importance. In addition, the analytical 
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expressions for the motions wovild also contain a constant term which 
represents the steady-state value of p and Ax if the system is stable. 
A good approximation to the time history could therefore be expressed in 
general form for po = -1.0, -I.5, -2.5, and -5.0 as 


or = Kq + K^e^'^sin(a5t + e) (5) 


For Pq = -1.86, -2.0, and -2.53 of case (b), and Pq = -2.0 of case (a), 
the expression is 


p(t) Ax(t) Alt Apt .. . 

— or — = Ko + Kpe i + K^e ‘ (4) 

Oq cxo ^ ^ ^ 

and for pq = -1.86 and -2.33 of case (a), 

PM or Mil = K.t + K.(e^^ - 1) (5) 

do OQ ^ 0 . 


Since the high-freqviency mode is omitted 


and ( 5) ^ the values of 


P(t) 


and 


Ax(t) 


OQ Oq 

condition of being zero at t = 0. 


from expressions (3)^ (4), 
will not satisfy the initial 


By differentiating expressions (3), (4), and (5) and setting them 
equal to zero, the maximm values of p and Ax and the time at which 
they occur could be easily determined. In these calculations, p and Ax 
occurred at about l/4 and l/2 of the period^ respectively. 


Approximation of the Amplitude Coefficients and Phase Angle 

The amplitude coefficients Kq, Kp, K^, Kj^, and K^, and 

phase angle e in expressions (5), (4), and (5) can be calculated directly 
by using the Heaviside expansion theorem given on page 45 of reference 6. 

A detailed examination of the exact expressions for the amplitude coef- 
ficients and phase angle indicated that many of the terms appearing in 
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the expressions had a negligible effect on the resialtant magnitudes. 
Thus the following simplified expressions were derived which result in 
good agreement with the exact values. 


In the expressions for 


P(t) 

°0 


Kq = — 
^ E 


MaNr 


% iz 


mV 


~ - Ix j 




Po 


(a^ + o^)(c - 2tu^) 




.0^ - 


Ma 


PO' 




Ma Nr 


Nr 


“ + 0^1 -2a + — + 

iy Iz \ 




. . „;!(. 2 a + ^ + 31 - rs 


€ = 


tan-l i 


Nr Mq\ La/Iy - Ix VlZ-^ 


mV' 


IPO' 


CD 


Iy 0 I Iz A Iy ; 




K2, and = 


%\iz ^ 


mV Ai 


2 CA + D 



f 

r- 


So , 

La 

Mq Nr ^ 3/% " ^x\ 4z " Lx\ 

Iy Iz ° V Iz / V Iy / 


D 

Iy niV 

4 
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In the 


where 


C2 


Cl 


K2, K5, 


expressions for 


a(t) 

°0 


Ko = - 


E 


% ~ ^ X] 


2 _ ipz - Ix' 

’ iz‘ 


% ^z 


Ki = \Jc2^ + Qi^ 


e 


tan-1 


Cl 



cuSCc - 2 oi2') 







and = - 



A(2CA + D) 



. %/^z - 

Iz\ % 



+ 5!a!k 

ly Iz 
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The C, D, and E appearing in the preceding expressions are the 
coefficients of the characteristic eqxiation. A compeirison of the motions 
obtained by using the exact and simplified coefficients is shown in fig- 
\ires 2 to 4 to be very good. 


Effect of Assuming pq Constant at t = 0 


In general, the rolling velocity builds up exponentially in response 
to an aileron input instead of being constant at t = 0 as assumed in 
this analysis. . In order to determine the effect of the assumption on 
the motions in p and Ax, an analog study was conducted assuming the 
rolling velocity reaches its steady-state value exponentially and the 
resultant motions in p and Ax were compared with, the motions in p 
and Ax obtained on the assumption that pq is constant at t = 0. 

Equations (la) to (id) were used in the analog study auid, in addition, a 


simplified roll equation 



was introduced. 


This equation 


states that p reaches its steady-state value, Pq, according to the 
expression p = Po(e^^ " l)* ^ general, the resxiltant motions obtained 


in the analog, study were very similar to the motions shown in figures 2 
and 5« The only significant difference noted was that different peak 


values of 


Ax 

aO 


and 


A. 

OO 


were encountered. 


Ax 


between the maximum values of — 

OO 


and 


OQ 


Table III presents a comparison 
obtained from the analog study 


and from figures 2 and 3, for pq = -1.0, - 1.5^ and -3»0. It is seen from 
the values presented in table III that higher maximum values occiar for 
Pq = -1.0 and -1.5 for the condition of constant Pq at t = 0 whereas 


the opposite is true for pq = '- 5*0. These resvilts can be explained by 

the fact that, during the first several seconds of the transient motion, 
the avergige vsilue of po in the analog study is smaller than under the 
assumption of constant Pq. Hence in the analog study, for Pq = - 1.0 

and -1.5 ^^he airplane may be considered to be located farther away 
from the stability boundary, but for pq = -3.0 the airplane is located 
closer to the stability boundary. 


Calculations of Motions for an Airplane Performing a 360° Roll 

The analysis presented thus far is applicable to the condition of the 
airplane performing a continuous rolling motion. Of particular interest 
are the motions of the airplane when the pilot performs a 360° roll. These 
motions can be approximated by first determing the values of 3, r, a, 
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and q at the time the airplane has rolled through 560°; that is, at 

Pif 

t = — . At this value of t, the constant rolling velocity is returned to 

Pq = 0 and the motion originally described by four degrees of freedom is 

now separable into its lateral and longitudinal values, the former repre- 
sented by the yawing and sideslipping equations and the latter represented 
by the pitching -moment and normal -force equations. Thus, the values of p, 

r, a, and q at t = — are the Initial conditions required to calculate 

PO 

the motions subsequent to t = — The values of p and a can be deter- 

Pq 

mined from equations (3), ( 4 ), or (5). From equations l(c) and 1 (d), the 
following equations are obtained: 


and 


, Yp 

r = — P + PoOO + Po ^ - P 


-■a 


q = — Al + PnP + a, 

mV ° 


( 6 ) 


(7) 


A good approximation of r and q at t — — is obtained by neglecting 

the terms p and a in equations (6) and (7). The following table com- 

pares the values of r and q at t = — for Pq = - 1 . 5 ^ -l- 7 > -5*0 

PO 

when p and a are Included in eqviations (6) and (7) and then deleted 
from equations (6) and ( 7 )* 


' Po 

P and a 

included 

P and a 

deleted 


q. 

r 

q. 

-1.5 

- 0.20 

0.20 

-0.19 

0.19 

-1.7 

C\J 

H 

1 

•35 

■ 

• 

1 

.32 

-3.0 

• 23 

.30 

.29 . 

.40 . 
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Effect of Gravity Terms on the p and 6a Motions 


The effect of the gravity components which would appear in the Y-force 
(eq. (ic)) and Z-force (eq. (id)) eqioations have been neglected in this 
analysis. In order to determine the effect' of the gravity terms on the p 
and Ax motions, an analog study was conducted using equations (ic) and 
(id) with the gravity terms included and then deleted from the equations. 
The airplane was assumed to be performing a 560 ° roll for rolling veloc- 
ities of pq = -I.5, -1*7j - 5 -b. A comparison of the maximum values 

of p and Ax obtained with and without the gravity terms is presented in 


table IV. The two values given for 

/Vy, 

respond to the maximum values of — 

OQ 


Ax ^1 Ax 2 

— , that is, — — and , cor- 

oq Oq Oq 

obtained during the initial part 


of the transient motion ^Pq = constant) and during the recovery part of 
the motion ^Pq = oj, respectively. The comparison shows that with the 

gravity terms included the maximum values of p are about 2 ° greater 
and the maxinium values of Ax increase by ^out 1° for ^ initial _oq _ 

of 5°. 


CONCLUDING RMARICS 


From the analysis presented in this paper it appears that the tran- 
sient motion in angles of attack and sideslip during a constant rolling 
maneuver consists chiefly of either an oscillatory mode or two aperiodic 
modes. Approximate expressions, are derived, for the determination of these 
modes as well as the modal coefficient corresponding to each mode . Inclu- 
sion of the derivatives CYo (side force due to sideslip) and Ct (lift- 

p a 

curve slope) increases the total damping of the system but does not 
markedly affect the transient motions. 

The sole input considered in this paper is the term PqOq (the 

product of rolling velocity and trim angle of attack) in the side-force 
equation. A similar analysis should be conducted for unit yawing moment 
and pitching moment inputs. With this information, the motion in angles 
of attack and sideslip due to pilot-applied rudder and elevator could be 
calculated. 


Langley Aeronautical Laboratory, 

National Ad'visory Committee for Aeronautics, 
Langley Field, Va., May 22, I 956 . 
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APPENDIX 

METHOD OF APPROXIMATING THE ROOTS OF THE CHARACTERISTIC EQUATION 

k ^ p 

The fourth-degree characteristic eqviation AA + + CX + DX + 

E = 0 may be factored into two quadratic equations ^X^ + aQ_X + 

^X^ + apX + = 0 where 

A = 1 

B = aj^ + ap 
C = bi + b2 + a]_a2 
D = apb]^ + a2^b2 
E - b^^bp 

An examination of the quadratic equations from which the exact roots 
presented in tables Il(a) and Il(b) were obtained indicated that bx » b2 
and that ax and ap are of the same order of magnitude and much smaller 
than bx* Thus, one may write 

A = 1 

B = ax + a2 

' C = bx 

, D = a2bx 

E = bxb2 
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The solution of these equations yields 




— - The -roots of- + a^^A + b]; = ~0 - approximate the high-frequency 
mode whereas the roots of A^ + a 2 A + b 2 = 0 approximate the remaining 

two roots of primary importance in the inotion calculations presented 
herein. A comparison between the exact roots and the approximate roots 
presented in tables Il(a) and Il(b) indicate that the agreement is 
excellent . 
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TABLE I 

MASS AND AERODYNAMIC CHARACTERISTICS OF AIRPLANE 


ni^ • • • • 

l x, slug-ft2 . . 

ly, . . 

sliig-ft^ . . 

Crip^ radian 

Cmq^ radian 

Cm^, per radian 
Cnp^ per radian 

Cyp> per radian 
Cl^, per radian 

Ct , per radian 
^P 

V, ft/sec . . 
sq ft . • . 

C y ft « • • • 

b y ft . • • • • 

q, Ib/sq ft 


745 

10,976 

57,100 

64,975 

-0.095 

- 3.5 

-0.36 

0.057 

-0.28 

3.85 

-0.255 

691 

377 

197 
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<11 


Figiire 1.- Boiondaries in the plane which define regions of 

aperiodic divergence for exaii 5 )le airplane . 



CONFIDENTIAL 


2k 


CONFIDENTIAL 


NACA EM L56FO4 



(a) Pq = -1.0. 


Figure 2.- Time histories of 


P/oo 


and Ax 


/«0 


for case a. 
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(b) Pq = -1-5- 


Figxore 2.- Continued. 
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Figure 2 .- Continued. 
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Figure 2.- Continued. 
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(e) Pq = -2.55. 
Figure 2 .- Continued. 


CONFIDENTIAL 






NACA RM L56f04 


CONFIDENTIAL 


29 



(f) Po = -2.5. 

Figure 2.- Continued. 
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(g) Pq = -5-0. 
Figiire 2.- Concluded. 
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(b) po = -1.5. 

Figure Continued. 
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(d) Pq - -2.0- 
Figure 3 *“ Continued. 
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(e) Pq = -2.53- 

Figure 3 .- Continued. 
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t , sec 


(f) Pq = 

Figure 5 *.- Continued. 
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(g) Pq = -5-0. 
Figure 5»- Concluded. 
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Figure 1<-.- Time histories of P/oq and for Cjj^ = -O.O9 and 

Cjip = O.llL and Pq = -I.O. 
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